Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 January 2010 



(MN M£X style file v2.2) 



Testing adiabatic contraction with SDSS elliptical galaxies 



A.E. Schulz 1 *, Rachel Mandelbaum 2 , Nikhil Padmanabhan 3 

1 Institute for Advanced Study, Einstein Drive, Princeton NJ 08540, USA 

2 Department of Astro-physical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA 
^Department of Physics, Yale University, New Haven, CT 06511 



5 January 2010 



o 

(N 
£ 

in 

O 

u 

Oh 

O 

> 
O 
VD 
(N 
(N 



o 

> 

i INTRODUCTION 



X 



ABSTRACT 

We study the profiles of 75 086 elliptical galaxies from the Sloan Digital Sky Survey 
(SDSS) at both large (70 — 700 hj , kpc) and small (~ 4 h^kpc) scales. Weak lensing 
observations in the outskirts of the halo are combined with measurements of the stel- 
lar velocity dispersion in the interior regions of the galaxy for stacked galaxy samples. 
The weak lensing measurements are well characterized by a Navarro, Frenk and White 
(NFW) profile. The dynamical mass measurements exceed the extrapolated NFW pro- 
file even after the estimated stellar masses are subtracted, providing evidence for the 
modification of the dark matter profile by the baryons. This excess mass is quanti- 
tatively consistent with the predictions of the adiabatic contraction (AC) hypothesis. 
Our finding suggests that the effects of AC during galaxy formation are stable to sub- 
sequent bombardment from major and minor mergers. We explore several theoretical 
and observational systematics and conclude that they cannot account for the inferred 
mass excess. The most significant source of systematic error is in the IMF, which would 
have to increase the stellar mass estimates by a factor of two relative to the Kroupa 
IMF to ful ly explain the mass ex cess without AC, causing tension with results from 
SAURON ( Cappellari et al.||2006 ) . We demonstrate a connection between the level of 
contraction of the dark matter halo profile and scatter in the size luminosity relation, 
which is a projection of the fundamental plane. Whether or not AC is the mechanism 
supplying the excess mass, models of galaxy formation and evolution must reconcile 
the observed halo masses from weak lensing with the comparatively large dynamical 
masses at the half light radii of the galaxies. 
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In hierarchical models of structure formation, matter orga- 
nizes itself into gravitationally bound halos that are formed 
through mergers of smaller halos, and by accretion along fila- 
mentary large scale structure. The properties of these halos 
have been studied extensively in N-body Cold Dark Mat- 
ter (CDM) simulations. Dark matter halos have been shown 
to exhibit a universal density profile, of which the most 
common parametrization is referred to as the NFW profile 
(Navarro et al. 111997 1 . More recent studies have determined 



that dark matter profiles deviate on average slightly from 
this two parameter description, and are better described by 
the Einasto profile, which introduces one extra parameter 
that mitigates the tendency for the measured concentration 
to spuriously depend on the innermost radius in the fit to 
the density profile ( |Navarro et al.|2004| |Merritt et al.|2005 
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Gao et al.|2008 1. The regularity in dark matter halos is one 
of the most robust predictions of hierarchical structure for- 
mation. 

Galaxies are believed to form inside dark matter halos 
via the cooling and condensation of baryons to their centers. 
Hydrodynamic simulations of hot gas in halos indicate that 
as this process occurs, baryonic physics affects the shape of 
the dark matter profile. As the baryons cool to the center, 
they begin to dominate the mass in the halo interior. The 
dark matter feels the change in the gravitational potential 
and is drawn into the center, steepening the central pro- 
file. The response of the dark matter to the cooling baryons 
has been modeled as the adiabatic contraction (AC) of a 
series of spherical shells under the assumptions of spherical 
symmetry, circular particle orbits, and conservation of an- 
gular momentum ( Blumenthal et aL]|1986 l. This treatment 
has been modified to account for the orbital eccentricity of 
particles in Gnedin et al. (20041, and has been shown to re- 



© 0000 RAS 



2 Schulz et al. 



produce the profiles of halos in hydrodynamical simulations 
at the 10-20 per cent level. 

Realistically, the AC process is synchronous with other 
more violent processes of merging and accretion. Older 
galaxies, in which the gas has been converted to stars long 
ago, continue to suffer dissipationless merging that has been 



argued to remove the influence of AC (Gao et al. 2004 Jo 



hansson et al.|20 09 ) . If a galaxy has had some portion of its 



baryonic content shock heated and removed to the warm- 
hot intergalactic medium where it is no longer seen, then 
potentially the contraction of the dark matter halo will be 
smaller, while the halo mass will remain the same. The hy- 



drodynamic simulations of Rudd et al. ( 2008J) suggest that 
the AC of the halo persists at a level that affects the mat- 
ter power spectrum at halo-sized scales. On the other hand, 



Naab et al. (20071 suggest that the galaxy formation his- 



tory, and specifically whether recent growth is due to star 
formation versus due to mergers and/or accretion, influences 
whether adiabatic contraction occurs or whether the dark 
matter halo in fact becomes less concentrated. Due to the 
many theoretical uncertainties, an observational handle on 
the mass distribution in both the inner and outer regions of 
the halo is crucial. 

The most reliable way to probe the dark matter is 
through its gravitational effects. Weak gravitational lensing 
(for a review, see Bartelmann fc Schneider||2001[ ) of back- 



ground galaxies by a foreground dark matter halo, or galaxy- 
galaxy lensing, is one robust way to constrain the halo den- 
sity profile. This method has been used to constrain halo 
density profiles outside of the regions where the influence of 
baryons is significant, for stacked galaxy samples (due to the 
lower signal) and for stacked and individual clusters (e.g. 
Dahle et al.|[2003l |Mandelbaum et al.||2006[ |2008a| |Okabe| 
et al.||20'09 Riemer-S0rensen et al.||2009 |. However, the in- 
ner regions of the halo where AC is important are typically 
inaccessible to weak lensing except perhaps for very deep 
space-based observations (e.g., Gavazzi et al.|[2007 : 



measure 



the weak lensing signal down to 10 h~ kpc). In this regime, 
the number of available background galaxies is small, confu- 
sion of light from the lens galaxy can make the measurement 
of background galaxy shapes unreliable, and the weak lens- 
ing approximation may fail entirely on small scales. Progress 
can be made by combining weak lensing with another probe 



such as strong lensing or kinematics (e.g., Kneib et al. 2003 



Gavazzi et al. 120071 [Limousin et al.|2007||Bradac et al.|200 8 
Corless et al.|2009||Oguri et al.|2009| |Umetsu et al .|2009| 



In this paper, we combine weak lensing observations 
in the halo, outside of the central regions (> 70 /i^kpc), 
with dynamical measurements that probe the central re- 
gions (< 10 ftf^kpc) of a stacked sample of elliptical galaxies 
from the Sloan Digital Sky Survey (SDSS). This approach 
allows us to constrain the density profile over a large dy- 
namic range of scales. An alternate approach is to combine 
strong and weak lensing measurements, as in |Gavazzi et al.] 
( 2007[ ); however, our approach allows the use of a statistical 
sample of elliptical galaxies, rather than being constrained 
only to those that are strong lenses. Finally, integral-field 
spectroscopy (e.g., de Zeeuw et al. 2002 Cappellari et al.] 



2006) can be used to probe the galaxy density velocity distri- 



bution and therefore infer the matter profile in detail within 
the effective radius, which in comparison with stellar pop- 
ulation synthesis models applied to the light distribution, 



can be used to infer the amount of dark matter within the 
effective radius. However, the connection to the profile on 
much larger scales is unclear, so it is impossible on the basis 
of IFU observations alone to distinguish between an adia- 
batically contracted halo versus a halo that is very massive 
overall. 

In section [2] we describe the sample of galaxies used 
in the analysis. In section [3] we begin with the dynamical 
mass measurement and then combine it with weak lensing 
observations. We further divide the sample along the mean 
size-luminosity relation and show that the scatter from this 
relation correlates with properties of the underlying dark 
matter profile. We conclude in section [4] Throughout, we 
use comoving distances unless specifically noted, with low- 
ercase r representing the three dimensional separation and 
uppercase R representing projected quantities. We have as- 
sumed a flat LCDM cosmology with Q. m — 0.3 and Qa = 0.7. 
Since we must compare stellar masses, which scale like h~ 2 
(for Hubble parameter Ho = 100/ikm/s/Mpc), with dynam- 
ical and halo masses, which scale like ft -1 , our results are 
by necessity dependent on the choice of h. We use h — 0.7 
for all quantities, and explicitly indicate their /i-dependence 
using /170; e.g., halo masses are indicated with K^qMq, stel- 
lar masses with K^§Mq, distances with /i^kpc, ratios of 
stellar to dynamical mass with h^g 1 , and similarly for other 
quantities. When relating results from other works that use 
h = 1, we simply use ft -1 kpc or /i _1 Mq. 



2 OBSERVATIONAL SAMPLE 

The lens and source samples that we use are from the SDSS 
(York et al. 20001, which has imaged roughly n steradi- 



ans of the sky, and followed up approximately one million 



of the detected objects spectroscopically jEtsenstcin et al 
200H |Richards et~al^|2002| |Strauss et al~1|2002[ ). The imag- 
ing was carried out by drift-scanning the sky in photomet- 
ric conditions (Hogg et al. 2001; Ivezic et al. 2004), in five 
bands (ugriz) (JFukugita et al. 1996; Sm ith et al.|2002' | using 



a specially-designed wide-field camera ( |Gunn et al. 



1998) 



All of the data were processed by completely automated 
pipelines that detect and measure photometric properties 



of objects, and astrometrically calibrate the data (Lupton 



et~aI1[200T| |Pier et al.|2003| |Tucke7eFa l. 2006; Padmanab 



han et al.|2008| ). T he SPSS was completed with its seventh 



data release, DR7 (Stoughton et al 
2003l|2004H2005l 



et al.||2006[ |2007 



2002 



Abazajian et al. 



Finkbeiner et al.|2004||Adelman-McC arthy 
Abazajian et al.||2009 l, from which 



2008 



we derive our sample, though it does not cover the entire 



2.1 Lens sample 

The lens galaxies are taken from the spectroscopic sample, 
for which we used the NYU Value Added Galaxy Catalog 
(VAGC, Blanton et al.|2005 |. For the portion of the analysis 
th at requires stellar masses , we have used the stellar masses 
of |Kauffmann et al.| ( |2003| ) updated for DR7Q 



.mpa-garching.mpg.de/SDSS/DR7/Data/stellarmass.html 
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The lens population is selected from the SDSS spec- 
troscopic sample using the following criteria. These crite- 
ria are designed to select elliptical or SO type galaxies with 
little or no rotational support that are the central galaxy 
in their dark matter halos, so that SDSS velocity disper- 
sion measurements can be reliably translated into dynamical 
masses, and the weak lensing signal more easily interpreted. 
All quantities are from the standard SDSS reductions (rerun 
137), unless specifically mentioned. 

• 0.02 < 2 < 0.35: We select lenses from the Main spec- 
troscopic sample with these spectroscopic redshifts. 

• g — r > 0.8 : This cut in rest-frame colour ensures that 
we are looking only at red galaxies with little to no active 
star formation. We have used model magnitudes fe-corrected 



to redshift 0.1 from the VAGC (Blanton & Roweis 20071 



The model magnitudes are extinction corrected using the 
reddening maps of Schlcgel ct al. (19981 and extinction-to- 
reddening ratios from Stoughton et al. (20021. 



• a v > 70 km/sec : This cut ensures that the stellar veloc- 
ity dispersion in the SDSS fiber can be reliably determined, 
given the resolution of the SDSS spectrograph. 

• Rgo/Rso > 2.6: This cut on the ratio of the Petrosian 
r-band 90 and 50 per cent light radii helps to ensure that 
we are selecting elliptical lens galaxies. The fact that the 
intensity profiles of ellipticals are more concentrated than 



those of spiral galaxies (Shimasaku et al. 20011 motivates 
this criterion. 

• b/a > 0.7: We select objects with large values of the mi- 
nor to major semi-axis ratio, obtained from model fits to a 
de Vaucouleurs profile. This cut will eliminate large edge-on 
spiral galaxies (which may pass our colour cuts due to inter- 
nal reddening), as well as elliptical galaxies with substantial 
rotation velocities that tend to flatten them. 

• Central: To ensure the lens galaxy is the central galaxy 
in its dark matter halo rather than a satellite galaxy, we as- 
sume that the brightest galaxy in a halo will be at the cen- 
ter. Thus we discard galaxies that have a brighter neighbor 
within a cylin der centered on the galaxy. Following |Reid fe| 
Spergel ( 2009 1, we select a cylinder of diameter 1.6 h~ Mpc 



(using h = 1 here for consistency with that work) and ± 864 
km/s along the line of sight. Due to fiber collisions in the 
SDSS, some of the galaxies that pass the cuts imposed on 
the Main spectroscopic sample lack redshifts. These galax- 
ies are not included in our lens catalog, but for the purposes 
of finding brighter neighbors, the redshift of the colliding 
galaxy is assigned to the object with no redshift to decide 
whether it falls into the cylinder. 

In addition to the selection made on the basis of lens galaxy 
properties, we also removed galaxies that lacked source 
galaxies in the background, that lacked stellar mass esti- 
mates, and that fell outside the mask of the random catalogs 
available in the VAGC. These cuts eliminated approximately 
15 per cent of the galaxies. 

For the bulk of our analysis, we have divided the 75 086 
lenses into three luminosity bins. The faintest bin contains 
half the lenses, the middle bin contains 1/3, and the bright- 
est bin contains 1/6 of the sample to roughly balance the 
signal-to-noise ratio (S/N) in the weak lensing measure- 
ments. Figure [T] shows the redshift distribution for each 
of these bins, and Fig. [2] shows the distribution of stellar 
masses. For each galaxy, we use the median of the PDF of 
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Figure 1. The redshift distribution of the the three luminosity 
bins used in this analysis. 
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Figure 2. The distribution of stellar masses. 



the stellar mass as our stellar mass estimate. A histogram 
of the comoving half light radii is plotted in Fig. [3] 



2.2 Source sample 

The source galaxies (with shape estimates) are taken from 
the SDSS photometric catalogs, with additional processing 



and selection cuts described in Mandelbaum et al. (2005a) 



This source sample contains over 30 million galaxies from the 
SDSS imaging data with r-band model magnitude brighter 
than 21.8, with shape measurements obtained using the 
REGLENS pipeline, including PSF correction done via re- 
Gaussianization (Hirata & Seljak||2003|), and with cuts on 



apparent size relative to the PSF designed to avoid various 
shear calibration biases. The overall calibration uncertainty 
due to all systematics was originally estimated to be eight 
per cent (Mandelbaum et al. 2005a I, though the redshift 
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Figure 3. A histogram of the (comoving) projected half light 
radii for each of the three luminosity samples. 



calibration component of this systematic error budget has 
recently been decreased due to the availability of more spec- 
troscopic data (Mandclb aum et al.|2008b| ). 

For the lensing analysis, a subsample of the total source 
catalog is used in probing the smaller radii, R < 100/i -1 kpc 
(described in section |3.2||. This procedure is a result of the 
analysis in |Mandelbaum et al. ( 2006[ l, which found that for 
bright elliptical lenses, use of the full source sample led to 
inclusion of satellites of those lenses that were intrinsically 
aligned towards their hosts, thus contaminating the lensing 
measurement. The LRG sources used at small separations 
were identified using colour cuts described in |Mandelbaum| 
et al. ( 2005a I to have very little contamination from galaxies 



below z = 0.4, thus minimizing contamination from intrin- 
sic alignments with our lenses, which are well below this 
redshift. 



3 RESULTS 

Table [l] summarizes all the properties of the lens sample 
used in this work. We now describe how we arrived at these 
numbers. 



3.1 Dynamical mass measurements 

We measure the dynamical mass at the 2D de Vaucouleurs 
scale radius for comparison with the weak lensing dark mat- 
ter profile measurement, described later in section |3~2| The 
dynamical mass is inferred by assuming a relation between 
the velocity dispersion of stars inside the portion of the 
galaxy covered by the SDSS fiber (3 arcsec in diameter), 
and the mass enclosed within the 2D de Vaucouleurs scale 
radius. The mass is related to the velocity of particles on 
circular orbits as 



v c (r P hys) = 



G M(r phys 

r phys 



(1) 



but relating the circular velocity to the radial velocity dis- 
persion a measured in SDSS requires a model for the stellar 



and dark matter mass distributions. The method for relating 
the two via the Jeans equation can be found in Pad manab-| 
han et al. (20041. We have assumed the stars follow a Hern- 
quist profile ( Hernquist|19 90 ) and the dark matter profile is 
an adiabatically contracted NFW profile, using the model of 
Gnedin et al. (20041. Since the sample is comprised of early 



type elliptical galaxies, we also assume that star formation 
has used all the available gas ( Read Sz Trentham|2005 Mor- 



ganti et al.||2006 Fukugita et al.||1998 l and we neglect any 
contribution to the mass from gas in the central region. We 
assume that there is no radial anisotropy /3 in the Jeans 
equation defined as 



1 



^=0 
vi 



(2) 



where v% is the radial velocity dispersion and v\ is the tan- 
gential velocity dispersion, which is not observable. Our re- 
sults are insensitive to this assumption, as discussed section 

El 

Fig. 3 compares the theoretical expectations for the ve- 
locity dispersion and circular velocity profiles for our mid- 
dle luminosity bin. The velocity dispersion profile changes 
by < 4 per cent from r/R^cv = 0.5 to 2.5, allowing us 
to extrapolate the velocity dispersion measurements from 
the fiber radius to the half-light radius. Furthermore, we 
find that the circular velocity is proportional to the velocity 
dispersion at the 5 per cent level, allowing us to robustly 
convert velocity dispersions into dynamical masses. At the 
projected half light radius, we find v c = 1.7a. Thus we can 
invert Eq. [T]to solve for the dynamical mass. 



A I, 



dy 



{l.laf R dc 

G(l+z) 



(3) 



Here we explicitly divide by 1 + z since our notation has 
RdeV in comoving coordinates. Fig. 2 of Padmanabha n et al.| 
(2004) shows that these curves are extremely insensitive to 



the value of the radial anisotropy parameter /3, even if /3 is 
not held constant with radius. The constant of proportion- 
ality is also insensitive (at ~ 5 per cent) to the stellar and 
DM mass profiles, and the details of the AC prescription, 



including no AC ( Padmanabhan et al.|2 004) . Variations in 
these model quantities cause M<j y to shift at the level of 
~ 10 per cent. We adopt this as the systematic uncertainty 
on the measurement of the total dynamical mass inside the 
SDSS fiber. 

Notice that the velocity dispersion a measured in SDSS 
is not measured at the projected half light radius Rdcv, but 
rather at the radius subtended by the SDSS fiber. We appeal 
to the flatness of u(r) demonstrated in Fig. [4] and assume 
that this measurement will equal the velocity dispersion at 
the half light radius. The angular diameter subtended by 
the SDSS fiber at the redshift of the middle luminosity bin 
is marked with a vertical line on the plot. 

We emphasize one detail: although the radius in Eq. [3] 
is the 2D projected half light radius, which is measured in 
SDSS by fitting to a de Vaucouleurs profile, the dynami- 
cal mass in Eq. [3] is an estimate of the 3D mass interior to 
that distance, not the mass in a 2D cylindrical projection 
through the halo out to that radius. Eq. [3] implies that at 
these radii the 3D dynamical mass is simply proportional 
to the radius (i.e., the profile is effectively isothermal), and 
we are choosing to make the measurement at the half-light 
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Table 1. Properties of the three lens samples studied in this work. The upper part of the table shows the result for the full sample, 
and the lower part shows the division of the sample along the mean size-luminosity relation (Smaller/Brighter or Larger/Dimmer than 
the mean relation). The notation elO is shorthand for (xlO 10 ). Luminosities are in absolute r-band magnitude, galactic extinction- and 
fc-corrected as described in section [2] The quantities M and c are fit to the weak lensing observations. The quantities z, M*, RdeV 
(comoving coordinates) and DM M^ y are lensing weighted averages over the lenses in the sample, and the Total M^ y is obtained by 
adding 0.42M* to the DM M<j y . M ntw is the mass expected from an NFW profile with mass M and concentration c. The dynamical 
masses and M n f w are measured inside i?dcV- The value h = 0.7 has been inserted into the distance modulus. The fraction of mass 
contributed by the stars at J?d c v is given by /*, and depends on the value of h. Errors on M, c, and M n f w are the 68 per cent CL 
obtained from bootstrap resampling of the NFW fits. Errors on M^ ot and i?dcV are measurement errors. Errors on the dynamical masses 
come from a 10 per cent estimate of systematic error in the mapping between the velocity dispersion and the mass interior to -RdcV- 




Figure 4. The circular velocity profile and radial velocity dis- 
persion multiplied by 1.7, which causes the two to agree near the 
scale radius. The parameters used to generate these curves are 
taken from the intermediate lens sample. Note that these are rel- 
atively flat functions over a wide range of radii. The vertical line 
marks the angular size of the SDSS fiber at the mean redshift of 
the middle luminosity sample. 



Figure 5. A histogram of the dark matter contribution to the 
dynamical mass for each of the three luminosity samples. 



the dynamical mass at the half light radius. A histogram of 
this quantity is plotted in Fig [5] 



radius in order to be able to estimate the stellar contribu- 
tion to the dynamical mass. The 3D stellar mass contained 
within the 2D projected half light radius i?d c v can be de- 
rived from the Hernquist profile to be roughly 42 per cent 
of the total stellar mass. We subtract this from the total 
dynamical mass to derive the dark matter contribution to 



3.2 Combining with weak lensing measurements 

Weak gravitational lensing probes the differential surface 
density of the dark matter associated with the lens galaxies. 
On scales probed here, this surface density can be attributed 
to the dark matter halo; on larger scales, due to large-scale 
structure. Thus, on the scales less than ~ 1 /i^Mpc, we can 



© 0000 RAS, MNRAS 000, 000-000 



6 Schulz et al. 



define the projected surface density via 
E(ii) = J p{R,z)Az. 



(4) 



Here p(r — y/R 2 + z 2 ) is the 3D density profile of the lens 
galaxies. 

Weak lensing is sensitive to the differential surface den- 
sity, defined as 



AE(7?) = E(< R) - £(i?) = 7tCR)E c , 



(5) 



(6) 



where the critical surface density E c is 

E ° 2 Ds 
c 4-kG D L D LS {1 + z L ) 2 

in comoving coordinates. Here, Dl and Ds are the angu- 
lar diameter distances to the lens and source galaxies, and 
Dls is the angular diameter distance between the lens and 
source galaxies. The right side of Eq. [5] shows how AE(i?) 
is related to the observable jt, the mean tangential shear 
of source galaxies behind the lens galaxies. This equation is 
true in the weak lensing limit (74 <C 1) in the presence of 
an axisymmetric mass distribution, which is achieved by our 
measurement procedure. 

The measurement is performed by stacking lens galaxies 
based on their luminosities. The stacking procedure allows 
us to achieve sufficient source number density to measure 
the shear signal, and therefore AE, of the composite ob- 
ject. To do so, we measure the tangential shear of all source 
galaxies in annular bins around the lens positions. The de- 
tails of the signal computation are described in appendix [A] 
We are implicitly assuming that on the scales used for our 
measurement, the shearing of background galaxies is caused 
only by the lens halos, and not by other nearby halos or 
infalling dark matter (in the language of the halo model, 
we neglect the two-halo term). For the few lenses that may 
be satellites despite our selection criteria designed to avoid 
them, we also neglect the contribution from the host halo. 
Furthermore, we neglect the radial window that enters the 
line-of-sight projection in Eq. [4] given that the window is 
far broader than the characteristic scale of our lenses. 

To generate model lensing signals that can be compared 
against the measurements, we fit the two-parameter NFW 
profile to the lensing data. The spherical NFW density pro- 
file is defined as 

Ps 



PnfwM 



(7) 



(r/r s ){l + r/r a ) 2 ' 

We can express it in terms of two parameters, the mass 
M within which the average density is equal to some over- 
density (here, 200 times overdense compared to the mean 
background density p), and the concentration, which is the 
ratio of the virial radius to the scale radius c = r v - 1T /r s . M 
and c can be related to p s and r s via 



M 



47T 'i 



47r / p(r) r dr 
'0 

3H 2 
l 8nG' 



P — Qmpc — ^r. 



(8) 



(9) 



Here Q m is the matter energy density today relative to the 
critical density, and H is the Hubble parameter. Our model 
signal can then be generated for a given value of M and c 
by projecting the profile along the line of sight (Eq. [4| and 
integrating to obtain AE (Eq. [SJ. 



Our use of a spherical NFW profile is justified because 
of the large number of lenses stacked, as verified by |Mandel-| 



baum et al. (2005b I. We use the NFW density profile rather 



than the Einasto profile, 



PEmasto(r) = p a e 



(-2/a)[(r/r,) -l] 



(10) 



but note that the two profiles are quite similar on weak lens- 
ing scales, though they are substantially different at smaller 
scales. One consequence of this choice is that the estimate 
of the concentration may potentially be systematically bi- 
ased in a way that depends on the innermost radius used in 
the fit. The impact of a systematic offset in concentration is 
discussed in section [3~4| However, as shown in Mandc lbaum| 



et al. (2008a I, fitting the galaxy cluster weak lensing profile 



on these scales for a concentration gave the same result to 
within several per cent when using an NFW profile and an 
Einasto profile. 

The fitting procedure involves two steps. In the first 
step, for each luminosity bin we fit the profile to the ob- 
served AE using the non-linear least-squares algorithm of 
Levenberg-Marquardt. Here, we use all data from 70 < R < 
2000 fe^kpc even though this may include scales on which 
lensing from large-scale structure is important. Our choice 
of minimum radius was motivated by the fact that ellipticity 
measurements of source galaxies closer than R < 70 /i^kpc 
to the lens center may be unreliable due to confusion of light 
from the lens, and smaller source counts. The resulting pro- 
file with (M, c) for each luminosity bin is used to estimate 
a virial radius for each bin. We then fit again, truncating 
the data at that virial radius to avoid significant contribu- 
tions to AE from material outside the lens halo (two-halo 
term) . This procedure results in fitting fewer radial bins for 
the faintest sample. The brightest and intermediate bins are 
fit in the interval 70 < R < 850 and the faintest bin in the 
interval 70 < R < 570 h^kpc. 

The lensing result and the NFW fit are shown in the 
top panel of Fig. [BJ The heavy solid line shows AE for the 
best fit NFW model, while the thin solid line shows the 
extrapolation of the NFW result to small radii. The points 
show the lensing measurements. The dashed lines shows a 
68 per cent confidence level estimate for the error on the 
NFW result for the middle sample. The errors are derived 
by bootstrap resampling 200 regions on the sky 1000 times, 
and recomputing the fit for each resampled dataset. The 
curves show the 16th and 84th percentile values of AE(J?) 
determined by rank-ordering the fits to the bootstrapped 
datasets. 

The bottom panel of Fig. [6] shows the (3D) integrated 
mass interior to the comoving radius r the NFW model of 
the middle sample shown in the top panel. For comparison, 
we also plot the enclosed mass for the Einasto profile, using 
the same mass and concentration parameters, and fixing the 



shape parameter following Gao et al. (20081. The horizontal 



scales in the two panels are the same. The total dynamical 
mass (upper point) and the dark matter contribution to the 
dynamical mass (lower point) are plotted at the position of 
the 2D de Vaucouleurs half light radius of the lens galaxies. 
The half light radius and dynamical masses shown are the 
lensing weighted mean value for the lenses in this luminosity 
bin. The lensing weights, which can be found in appendix 
|A"} are assigned to each lens-source pair, and include red- 
shift information and error on the shape measurement of 
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Figure 6. Top: Points are the weak lensing measurement. Solid 
lines show the differential surface density for the best fit NFW 
model. Dotted lines show 16 and 84 per cent contours from fits 
to bootstrap resampled datasets. Bottom: Solid lines show the 
integrated mass interior to r, and dotted lines show a bootstrap 
estimate for error in NFW. The dashed curve shows the result 
for an Einasto profile for comparison. The upper point shows 
the total dynamical mass, while the lower point shows the dark 
matter contribution to the dynamical mass after subtracting the 
estimated contribution from stellar mass. 



the source. To deduce the dark matter contribution to the 
dynamical mass, we subtract 42 per cent of the total stel- 
lar mass (the fraction of stars enclosed within a 3D radius 
corresponding to the 2D half light radius) from the total dy- 
namical mass. The error bar on 7?d c v comes from the mea- 
surement error. The error bar on the total dynamical mass is 
estimated to be ~ 10 per cent, from variations in the input 
modeling from Eqn. [3] and from differences between Rd e y 
and the portion of the galaxy covered by the SDSS fiber. 
The error bar on the dark matter dynamical mass comes 
from combining this systematic error in the total dynamical 
mass with measurement errors in the stellar masses. 

Although the NFW curve has been extrapolated inward 
by more than an order of magnitude in transverse separa- 
tion, it is noteworthy that the total dynamical mass in the 
interior of the galaxy exceeds the NFW prediction by eleven 
standard deviations. Even when the contribution of stars 
has been accounted for, there appears to be over a factor 
of two more dark matter mass in the central region of the 



galaxy (four standard deviations) than what is seen in halos 
from N-body simulations of large scale structure. 

This result is consistent with earlier attempts to observe 
the density profile of elliptical galaxies. |Jiang fc Kochanek| 
(20071 used stellar dynamical and strong lensing measure- 



ments of 22 early-type galaxies from the Sloan Lens ACS 
Survey (SLACS) to constrain the total matter profile near 
the effective radius. These data alone are not enough to dis- 
tinguish between a model with and a model without adia- 
batic contraction; however, when the y include weak lensin g 
constraints on 100/i - kpc scales from Gavazzi et al. (2007), 



they find that the models with adiabatic contraction are 
strongly favored. For their analysis, they used the |Blumen-| 
thai et al. (1986) model for adiabatic contraction, which pre- 
dicts a more significant effect than the |Gnedin et al.| ( |2004[ , 
model used in our work. 



3.3 Comparison with adiabatic contraction model 

It is well known that baryonic physics is expected to affect 
the form of the dark matter profile in galaxy halos, so it is 
not entirely surprising that the profile departs from N-body 
predictions in a region where the mass is largely made up of 
stars. Exactly how the baryons affect the dark matter dis- 
tribution, and whether their effect is transient or enduring 
is however a matter of some theoretical debate. On the one 
hand, the dark matter is expected to fall into the cuspy deep 
gravitational potentials that are created when baryons cool 
and condense to form galaxies, a process known as adiabatic 
contraction (AC) that further deepens the central potential. 
In the simulations of |Rudd et al.| ( |2008| ), these steeper po- 
tentials are seen to persist, and the impact of the baryons is 
evident on the matter power spectrum, principally through 
the change in the dark matter profiles inside the virial radii 
of halos, particularly around cluster masses of W^h^Mo. 
On the other hand, there is some debate as to whether spiky 
central cusps will survive the onslaught of subsequent major 
and minor merger activity expected in hierarchical models 
of structure formation. In the simulations of Johansso n et al.l 
( |2009[ ), adiabatic contraction is seen initially, but is trans- 
formed to a cored central potential through gravitational 
heating from mergers. To date, there have been very few 
observational studies that can probe the dark matter halo 
over the large range of scales required to test the AC hy- 
pothesis. We make an attempt here by comparing the AC 
model of Gnedin et al. (2004) to the dynamical mass data, 
shown in Fig. [7] 

Each of the panels in Fig. [7] is analogous to the bottom 
panel of Fig. [6] but each quantity is now divided by the best 
fit NFW relation (solid curve in Fig. [6| so as to present the 
results on a linear scale. Note that the best fit NFW is dif- 
ferent for each panel. The 3d stellar fractions of the total 
mass at the half light radius are /» = 31, 35 and 38fef per 
cent (from bright to faint). The curves in Fig. [7] show the 



AC models of Gnedin et al. ( 2004 \ (dashed) and |Blumen 



thai et al. (19861 (dot-dashed). These are computed using 
the publicly available contra code ( Gnedin et al.|2004 i. To 
model AC, we use the mass and concentration from the weak 
lensing measurements, and make a few reasonable assump- 
tions. We assume that for these elliptical lens galaxies all 



of the gas has been converted into stars ( Read & Trentham 



2005 Morganti et al. 2006 Fukugita et al. 19981, so that 
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Figure 7. All curves and points are defined in the same con- 
vention as in the bottom panel of Fig. [fj] with the addition of 
a dashed line showing the predictions including AC. The non- 
AC NFW model has been divided out (solid line at 1) so as to 
present the results on a linear scale. The (top, middle, bottom) 
panel shows the (brightest, intermediate, faintest) lens sample. 



the baryon fraction is given by the stellar mass divided by 
the best-fitting dark matter halo mass. This assumption is 
robust for the fainter bins but may begin to break down for 
the brightest bin, for which the mass approaches the group 
scale. We also assume the baryons ultimately arrange them- 
selves into a Sersic profile with n = 4, i.e. a de Vaucouleurs 
profile, with scale radius J?d c v converted as in the appendix 



of Padmanabhan et al. (20041. 



It is remarkable that in all three luminosity bins, the 
model NFW profile with adiabatic contraction agrees with 
the inferred contribution of dark matter to the total dynam- 
ical mass. This AC model leads to enhancements of the dark 
matter contribution to the dynamical mass that are approx- 
imately factors of 4, 3.5, and 2 for the faint, middle, and 
bright samples respectively. The disagreement of the data 
with the extrapolated pure NFW profile (i.e., without AC) 
is significant at the 3.5a, 4<r, and 2.5a level, respectively. 
The apparent agreement between the total dynamical mass 
and the Blumenthal et al. ( 1986 1 model of AC is a coinci- 



dence; for this to be meaningful, none of the mass could be 
in stars. 



3.4 Potential systematic errors 

In this section, we discuss a number of systematic uncer- 
tainties present in this measurement that make the very 
close agreement between the AC model and the measure- 
ment mildly surprising. 



3.4-1 Initial mass function 

The most significant concern for potential systematic er- 
ror lies in the determination of the stellar mass. The stellar 
masses are computed assuming a Kroupa Initial Mass Func- 
tion (IMF, Kroupa|2001 1, which is calibrated by measuring 
the mass function of stars in the solar neighborhood. While 
it is expected that some level of churning will bring stars 
from orbits in all parts of the galaxy close to the solar neigh- 
borhood, the degree of churning is very uncertain. More im- 
portantly, galaxies of different ages and types have different 
metallicities, and the shape of the IMF may well depend on 



metallicity. In Kauffmann et al. ( 2003 1, the sensitivity of the 



stellar masses to the choice of IMF was investigated, and it 



was shown that assuming a Salpeter ( Salpeter 1955 1 IMF 
with a lower mass cutoff at O.IM© systematically increased 
the stellar masses by approximately a factor of 2. However, 
the amount of increase in the stellar masses will quantita- 
tively depend on how the low mass end of the (divergent) 
Salpeter IMF is regularized. It is interesting to contemplate 
that aside from the question of AC, this measurement of the 
dynamical mass can be thought of as a constraint on the 
IMF, in the sense that the stellar mass cannot exceed the 
total dynamical mass. 

To fully explain the observed dynamical masses with no 
adiabatic contraction of the dark matter, the stellar masses 
would have to increase by a factor of 1.9, 2.3 and 2.2 for the 
bright, middle and faint samples. However, |Cappcllari et aL] 
(2006) present a challenge to this alternative explanation. 
They use detailed integral field spectroscopy observations 
to infer dynamical masses within the effective radius for 25 
nearby E/S0 galaxies, and compare them with the predicted 
stellar masses from stellar synthesis models. They find that 
for some of the galaxies, the predicted stellar masses with a 
Salpeter IMF exceed the dynamical masses, which suggests 
that if we require the same IMF for each galaxy, then the 
Salpeter IMF is too bottom-heavy. In contrast, the Kroupa 
IMF gives stellar masses that are 30 per cent lower than with 
the Salpeter IMF (for their choice of cutoff mass), causing 
all the stellar masses to be less than or equal to the inferred 
dynamical masses. Thus, the conclusion that the IMF is very 
bottom-heavy (instead of accepting the AC hypothesis) may 
be difficult to reconcile with these IFU observations of ellip- 
tical and SO galaxies. 



3-4-2 Other stellar mass uncertainties 

There are other systematic effects in the stellar mass es- 
timates that are expected to be subdominant to the IMF 
uncertainties. For example, the effects of aperture bias (the 
fact that the M* / L are derived from the region of the galaxy 
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covered by the fiber, not the whole galaxy) are estimated to 
be smaller than that of IMF uncertainties. Also, aperture 
bias is more problematic for galaxies with both differently- 
coloured bulge and disk contributions than for the pure el- 
lipticals of which our sample is composed. Sec Kauffmann 
et al. (20031 for a more full discussion of all stellar mass- 



related systematics. 



3.4.3 Initial conditions 

Another systematic uncertainty relates to the conditions un- 
der which adiabatic contraction initially occurred. The pro- 
cess of AC happens while the baryonic component is still in 
gaseous form, whereas the measurement is being made on 
galaxies that have converted nearly all of the gas into stars. 
Presumably several (major or minor) mergers have occurred 
in the interim, whose qualitative effect will be to raise the 
mass and lower the concentration of the halo. If the baryon 
fraction of the merging objects differs from that of the par- 
ent halo, the baryon fraction will also have evolved since 
the time AC was occurring. To quantify how sensitive the 
AC result is to changes in the input parameters, we show 
in Fig. [8] the effects of a 30 per cent variation in halo mass, 
concentration, scale radius of the baryon component, and 
baryon fraction. 

The CONTRA code accepts the baryon mass fraction 
as an input, and returns the pre-contracted and post- 
contracted curves in units where the mass has been scaled 
to the NFW virial mass, which we measure in the lensing 
analysis. We define the baryon fraction as = M,/M (as- 
suming that all the gas has been converted to stars Q The 
top two panels of Fig. 8 show the result of varying the halo 
mass M, in one case fixing the baryon fraction (top) and 
in the other case allowing the variation of M to affect the 
baryon fraction (second panel). In the top panel, changing 
the mass simply scales the curves up and down. The second 
panel shows that the change from scaling the mass in the 
top panel is nearly exactly offset at these radii by the effect 
of decreasing the baryon fraction. This result can be under- 
stood by comparing the top panel with the bottom panel 
(where only the stellar mass is scaled). A decrease in stellar 
mass implies a decrease in the baryon fraction, and since the 
deviation from the fiducial model is roughly the same size 
in the top and bottom panels, the effects cancel one another 
out in the second panel. One other noteworthy feature in 
Fig. 8 is that the effect of increasing the baryon scale radius 
changes the level of AC in the opposite sense than the other 
4 curves. We conclude that there is a modest sensitivity in 
the model to these quantities that could change the predic- 
tion at a level comparable to the measurement error of the 
dynamical mass (indeed this is the reason we selected a 30 
per cent variation). This sensitivity is not so large as to in- 
validate the result for reasonable (tens of per cent) shifts in 
these parameters. 



2 Please note this is distinct from /„ in Tabic [7J which indicates 
the fraction of the total dynamical mass that is in stars within a 
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Figure 8. Each panel shows AC of the best fit NFW model (Solid 
curve), the total dynamical mass (upper grey point) and the dark 
matter contribution to the dynamical mass (lower black point). 
The dashed [dot-dashed] curve shows the impact on AC of in- 
creasing [decreasing] the halo mass (upper), concentration (upper 
middle), scale radius of the baryons (lower middle), and stellar 
mass (lower). All results have been divided by the non-AC best 
fit NFW model. 



3.4-4 Dynamical mass assumptions 

It is quite probable that the lens sample violates the assump- 
tion that the radial anisotropy parameter ft is zero. It has 



already been shown in Padmanabhan et al. (2004 1 that the 



measured velocity dispersion is insensitive to the value of the 
radial anisotropy parameter. However it is also worth not- 
ing that for this type of galaxy, the anisotropy is much more 



likely to be radial rather than tangential ( |Cappellari et al. 
2007 1 . The impact of assuming no radial anisotropy for such 



sphere of radius equal to the half light radius R, 



dcV- 



objects will be to systematically underestimate the dynami- 
cal mass within the half light radius. A similar systematic ef- 
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feet that causes the dynamical mass to be underestimated is 
if the galaxies have some moderate rotational support. Both 
of these effects strengthen the conclusion that there appears 
to be more mass than predicted by simple extrapolation of 
the NFW profile; our dynamical mass measurements can be 
considered a lower bound. 



3.4-5 Spread in the properties of the lens subsamples 

Finally, our analysis uses relatively broad luminosity bins 
(because of the need for large lens samples to measure the 
lensing signal at reasonable S/N) . Unfortunately, this means 
that a large amount of averaging takes place, which could 
muddy the interpretation of the results. We have attempted 
to ameliorate this effect by computing lensing-weighted aver- 
age quantities (stellar masses, dynamical masses, half-light 
radii, and so on) to compare against the best-fitting dark 
matter halo mass. This procedure should help make the dif- 
ferent kinds of measurements more comparable. However, 



Mandelbaum et al. (2005b I show that for lens samples with 



halo mass distributions that are more than a factor of five 
wide, the best-fitting halo masses lie somewhere above the 
median and below the mean, an effect that can be as large as 
several tens of per cent in extreme cases. Since Fig. [8] shows 
that changing the halo mass and other properties by 30 per 
cent does not change the AC predictions by more than about 
la, we anticipate that our neglect of the bin width (which is 
unknown, in the case of the halo mass) will not substantially 
alter our conclusion that the AC model is consistent with 
our findings. 



3.5 Splitting along the size-luminosity relation 

It is interesting to ask if the level of excess mass in the galaxy 
interior correlates with other galaxy properties. There is ev- 
idence that a galaxy's environment and merger history di- 
rectly impact the concentration and relative distributions of 
baryons and dark matter ( Abadi et al.||2009 and references 
therein). It has also been suggested that multiple major or 
minor mergers may erase some of the effects of adiabatic 
contraction ( Johansson et al.|2009 l. Since elliptical galaxies 
lie tightly on a fundamental plane, it would also be interest- 
ing if small departures correlate strongly with the form of 
the galaxy's dark matter profile. Fortunately, there is enough 
signal in the data to support one further level of subdivision 
in the analysis. We have opted to examine these questions by 
dividing each luminosity bin along the mean size-luminosity 
relation, which is a projection of the fundamental plane. The 
mean relation in our sample is 



log(-Rdev) = 0.722 + 0.598 * ((-20.44 - M r ) /2.5) 



(11) 



where RdcV is in units of h? kpc. The results are presented 
in Tableland Fig. [9] The upper set of curves (dashed lines) 
show the departure from NFW (solid line at 1) due to AC for 
the smaller/brighter sample. The lower set of curves (dot- 
dashed lines) show AC for the larger/dimmer sample. The 
points on the left (right) correspond to the smaller/brighter 
(larger/dimmer) sample. Table [I] indicates that for the faint 
sample, we are unable to simultaneously constrain M and c 
for these subsamples divided along the mean size-luminosity 
relation. Therefore, in the bottom panel of Fig. |9j we adopt 
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Figure 9. Three luminosity bins are shown from brightest (top) 
to dimmest (bottom panel). In each panel, curves on the top (bot- 
tom) show the AC results for lens subsamples that are smaller and 
brighter (larger and dimmer) the mean size luminosity relation. 
The solid lines with dotted error estimates show the enclosed mass 
for the NFW model that best fits the weak lensing data. The to- 
tal observed dynamical mass (upper grey point) and the dark 
matter contribution to the dynamical mass (lower black point) 
plotted at the position of the half light radius. Points on the left 
(right) are the smaller/brighter (larger/dimmer) sample. Because 
the mass and concentration cannot be independently constrained 
for the two subsets of the faintest sample, M and c are fixed to 
the values of the composite faint sample (top level of Table [TJ . 



the M and c from the composite sample (row 3 of Table [T]) . 
As a result, for the top two panels, the AC curves and data 
points for the two subsamples have been divided out by a 
different non-contracted NFW model; on the bottom panel, 
by the same model. 

Table [I] shows that the halo masses and half-light radii 
of the smaller/brighter (S/B) sample are smaller than those 
of the larger/dimmer (L/D) sample, but the stellar masses 
and concentrations of the two populations are nearly the 
same. Fig. [9] indicates that the total dynamical mass at the 
half light radius is much more enhanced above the NFW 
prediction for the S/B sample than for the L/D sample. 
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There are competing effects driving this trend; it is helpful 
to refer to Fig. [8] in deciphering the result. Since the stellar 
masses are equivalent, the smaller halo masses for the S/B 
sample imply a larger baryon fraction ft, than in the L/D 
sample. Just as in panel two of Fig. |SJ the increase in mass 
is offset by the decrease in the baryon fraction in compar- 
ing the S/B to the L/D sample. The concentrations of the 
two samples do not change very much, but the significant 
change in the baryon scale radius -RdeV drives the overall 
trend. Since the S/B sample has a larger baryon fraction, 
and more of the baryons are concentrated toward the center 
of the halo, the effects of AC are more pronounced than in 
the L/D sample, even though the L/D sample has higher 
overall halo mass. There is also a higher stellar mass frac- 
tion /„ = (0.42M„)/Mj y °* at the half light radius for the 
S/B sample than for the L/D, although part of this effect 
could be accounted for by non-homologous profiles (P. v. 
Dokkum, priv comm). 

An alternate interpretation of this result might be that 
the stellar IMF for the S/B sample has more mass in small 
faint stars than the IMF for the L/D sample (see the discus- 
sion of systematic error in Section 3.2 above), or potentially 
some combination of these two effects. The IMF interpreta- 
tion requires us to discard the idea of a universal IMF for 
all elliptical galaxies. 

Since it is likely that the S/B sample had its stellar 
population in place at earlier times, and that the L/D has 
potentially experienced merger activity more often or more 
recently, it would be very interesting in the future to corre- 
late the level of AC with other indicators of age or merger 
activity. It is also interesting to note that the L/D subsam- 
ples within each luminosity bin have consistently larger halo 
masses than the S /B subsamples, despite the fact that their 
stellar masses are very similar. This result suggests that dif- 
ferent star formation histories (that led to the different ap- 
pearance of the stellar components at fixed stellar mass) 
also correlate with different accretion histories for the dark 
matter halos. 



3.6 Comparison with previous work 

Here we discuss how our results compare with previous, re- 
lated observations. We first compare with Padmanabhan 
et al. (20041, which is relevant due to the similar sample 



definition and methodology. Since that work did not have 
any observational method of constraining the profile on large 
scales, we simply compare against their determination of the 
3d stellar mass fraction /« within the half-light radius. The 
relevant comparison is against figure 8 of that paper, which 
shows Mdy/M* (i.e., our /<T ) as a function of the galaxy ap- 
parent size, for galaxy subsamples split by stellar mass. For 
our three relatively broad luminosity subsamples, we have 
typical values of half-light radius range from 4 to lS/i^kpc 
and typical 1//* ranging from 2.6 to 3.2/170; Padmanabhan 



et al. (20041 find values of 1//, that range from 2-4 (de- 



pending on the stellar mass, in narrower bins than ours), 
which is consistent with our findings. 



The comparison withlGavazzi et al. (20071 is less triv- 



ial due to the very different methodology. In that paper, 
they model the strong and stacked weak lensing signals for 
SLACS lenses jointly using the sum of an uncontracted NFW 
and a stellar Hernquist profile with a free M*/L (without 



any stellar mass estimates from stellar population synthe- 
sis methods). With this procedure, they quote a 3d stellar 
mass fraction within the effective radius of /« = 0.73 (inde- 
pendent of h because of how the modeling is done). However, 
|Jiang fc Kochanek (20071 model the same SLACS sample, 
and find that if they allow for the possibility of adiabatic 
contraction, the resulting M*/L is 70 per cent of the value 
when adiabatic contraction is ignored (see their figure 6). 
This decrease results from a degeneracy in the fits: an un- 
contracted DM profile requires more mass to be attributed 
to stars than a contracted DM profile, in order to account 
for the total mass. Furthermore, they find that the fit with 
AC (and lower value of M*/L) is preferred once they include 
the weak lensing in the modeling. Consequently, the |Jiang fc| 
Kochanek] (|2007[) res ults suggest that this value of /, from 
Gavazzi et al. (20071 should be reduced to ~ 0.5, which is 



higher than our values of 0.31-0.38, but not very signifi- 
cantly once the different sample selections (strong lensing 
versus optical) and redshift ranges are taken into account. 
Finally, we compare against the results from Cappcllari 



et al.| ( |2006| ), who find that for 25 SAURON E/S0 galaxies, 
the assumption of a Kroupa IMF yields /', ~ 0.7 (median 
value, independent of h). While this number is significantly 
higher than ours, the selection criteria and redshift range 
of the galaxies are significantly different, so it is difficult to 
derive any conclusions from this comparison. 



4 CONCLUSIONS 

We have combined weak lensing and velocity dispersion ob- 
servations to study the dark matter profiles of elliptical 
galaxies from the SDSS. The radial profile of the weak lens- 
ing shear is consistent with that expected from an NFW 
halo profile. We have fit two NFW parameters to the weak 
lensing data (mass and concentration), and extrapolated the 
profile inward to smaller radii not accessible to weak grav- 
itational lensing. We have deduced the dynamical mass at 
these smaller radii by measuring the velocity dispersion of 
the stars within the half light radius of the galaxy. We com- 
pare the measured dynamical mass to the extrapolation of 
the NFW profile, and find that there is a significant excess 
of mass in the interior. Using estimates of the stellar mass 
of the galaxy, and assuming that all non-stellar mass is dark 
matter, we find that the dark matter contribution to the 
dynamical mass is still far in excess of the NFW prediction. 
This result is in support of the model of adiabatic contrac- 
tion (AC), which predicts that the cooling and condensation 
of baryonic material will deepen the gravitational potential 
well of the galaxy, and pull the dark matter towards the 
halo center. These results suggest that the effects of AC are 
stable to the subsequent bombardment of major and minor 
mergers suffered by these objects since the time their gas 
cooled. 

We compare the observation to a theoretical model of 



adiabatic contraction (Gnedin et al. 20041 and find good 



agreement between the two. However, we acknowledge that 
systematic uncertainties in the determination of the stellar 
mass make it difficult to prove that the excess mass at the 
half light radius is indeed dark matter. An alternative inter- 
pretation may be that the excess mass is comprised of some 
other dark baryonic form, for example it may indicate an 
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IMF with more of the mass concentrated in dim, low mass 
stars relative to the Kroupa IMF used in our stellar mass 
estimates. In order to fully explain the observation without 
adiabatic contraction, the stellar masses would need to in- 
crease by a factor of two. A change of this magnitude in the 
stellar masses would be difficult to reconcile with recent re- 
sults from the SAURON collaboration. They conclude that 
a 30 per cent increase in stellar mass (from masses obtained 
with the Kroupa IMF) is inconsistent with their dynamical 
data derived from integral-field spectroscopy. 

Finally, we divide the sample of galaxies along the mean 
size-luminosity relation and show that the scatter from the 
mean relation is related to the underlying properties of the 
dark matter profile. We find that due to a higher baryon frac- 
tion with more of the baryons concentrated toward the cen- 
ter, the smaller/brighter subsample experiences more lasting 
effects of adiabatic contraction than the larger/dimmer sub- 
sample. If adiabatic contraction is not the explanation, the 
most likely alternative is different IMFs for the two subsam- 
ples. Also, we find that the two subsamples have comparable 
stellar masses but different halo masses, suggesting that stel- 
lar mass does not trace halo mass, possibly due to different 
formation histories that result in differing dark and stellar 
components (with larger halo mass being associated with a 
dimmer, less concentrated stellar component). 

Ultimately, regardless of whether the excess mass is 
dark matter or baryonic, this set of observations indicates 
that the total mass at the half light radius far exceeds what 
would naively be expected from an NFW or Einasto pro- 
file. Contemporary models of galaxy formation and evolu- 
tion must accommodate this relatively high ratio of dynam- 
ical mass in the interior to total halo mass in early type 
elliptical galaxies. 
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APPENDIX A: COMPUTATION OF THE 
WEAK LENSING SIGNAL 

We compute the lensing signal in logarithmic radial bins 
from 30 h^kpc to 2.9 fe^Mpc as a weighted summation 
over lens-source pairs in the bin, using the following estima- 
tor: 



AE(R) 



E(ls) vri(ls) 
Is Visit 



(Al) 



where the factor of 27?. arises due to the response of our 
cllipticity definition to a shear. 1Z, known as the "shear re- 



0.87 (Bernstein & 



sponsivity," is approximately 1 — e^ n 
Jarvis|2002| ). 

The weight factors wi s assigned to each pair include 
redshift information and the error on the source shape mea- 
surement via 

E" 2 



Wis = 



+ o\ 



(A2) 



where o"1n, the intrinsic shape noise, was determined as a 
function of magnitude in [Mandelbaum et al. (2005a|, fig. 
3. The factor of E^T 2 downweights pairs that are close in 
redshift, so that we are weighting by the inverse variance of 
AE. 

The critical surface density, E c , is computed in two dif- 
ferent ways, depending on the source sample. For sources 
with r < 21, we use photometric redshifts, and require 
that the source photometric redshift be greater than the 
lens redshift plus 0.1. The procedure for correcting these 
E c estimates to account for photometric redshift error is 



described in Mandelbaum et al. (2005a 2008b I. For sources 
with r > 21, we use a redshift distribution for the sources, so 
that for a given lens redshift a single E c value is used for all 
sources at r > 21. Finally, for the high-redshift LRG source 
sample used only at transverse separations < 140/if 1 kpc, we 
use photometric redshifts to assign E c for each lens-source 
pair. 

There are several additional procedures that must be 



done when computing the signal (for more detail, see Man- 
delbaum et al. 2005a I. First, the signal computed around 



random points must be subtracted from the signal around 
real lenses to eliminate contributions from systematic shear. 
In practice, this correction is negligible on the scales used 
for this measurement. Second, the signal must be boosted, 
i.e. multiplied by B(R) = n(R)/n ran d(R), the ratio of the 
number density of sources around real lenses relative to their 
number density around random points, in order to account 
for dilution by sources that are physically associated with 
lenses, and therefore not lensed. 

In order to determine errors on the lensing signal, we 
divide the survey area into 200 bootstrap subregions, and 
generate 1000 bootstrap-resampled datasets. 



© 0000 RAS, MNRAS 000, 000-000 



